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Abstract 

We study the generating function of rooted and unrooted hyperforests in 
a general complete hypergraph with n vertices by using a novel Grassmann 
representation of their generating functions. We show that this new approach 
encodes the known results about the exponential generating functions for the 
different number of vertices. We consider also some applications as counting 
hyperforests in the fc-uniform complete hypergraph and the one complete in 
hyperedges of all dimensions. Some general feature of the asymptotic regimes 
for large number of connected components is discussed. 
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1 Introduction 

In this paper we shall be concerned mainly with the problem of evaluating the 
weight of rooted and unrooted hyperforests in the complete hypergraph with n ver- 
tices K,n when the weight of a hyperedge depends only on its cardinality. These 
questions are usually analyzed by using the exponential generating function and the 
Lagrange inversion formula [1,2], eventhough it seems that they have been posed and 
solved in the context of statistical mechanics [3]. But, at least in the case of ordinary 
graphs, the entropy of trees and rooted forests of a generic graph can be evaluated by 
using Kirchhoff 's matrix-tree theorem. For the case of unrooted forests a solution can 
be obtained by the use of a novel generalization of the Kirchhoff 's theorem [4], where 
the generating function of spanning forests in a graph, which arises as the g — > 
limit of the partition function of the g-state Potts model [5-8] , can be represented as 
a Grassmann integral involving a quadratic (Gaussian) term together with a special 
nearest-neighbor four-fermion interaction. Furthermore, this fermionic model pos- 
sesses an osp(l|2) supersymmetry. By applying this method the classical result [9,10] 
that the number of unrooted forests on the complete graph with n vertices for large 
n behaves asymptotically as n"~^-v/e can be recovered [11]. But also more detailed 
informations. For example in [12] the renormalization flow for unrooted forests on 
the triangular lattice has been analyzed. 

A further generalization has been achieved in [13], where, given a hypergraph 
G = {V, E) (that is, E is an arbitrary collection of subsets of V, each of cardinality 
> 2), by exploiting the underlying osp(l|2) supersymmetry, a class of Grassmann 
integrals permits an expansion in terms of spanning hyperforests. More precisely, let 
us introduce, at each vertex i G a pair of Grassmann variables ipi^ ipi, which obey 
the usual rules for Grassmann integration [14, 15]. For each subset A C V wc define 
the monomial ta = YiieA'^'i'^i^ each number A (in M or C), we define the 

Grassmann element 



and introduce a notation for the integral on all the Grassmann fields on the vertices 



Given arbitrary hyperedge weights {wA}AeE, the general Grassmann integral 
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( "partition function" ) 
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has a combinatorial interpretation in terms of spanning hyperforests. 

Special cases provide the generating functions for rooted and unrooted spanning 
(hyper) forests and spanning (hyper)trees. The generating function of unrooted span- 
ning hyperforests, with a weight wa for each hyperedge A and a weight A for each 
connected component is given by 



/ 



Vv{ip,ip) exp 



A ^ ipiipi + ^ wa/a 



(A) 
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Wa^ 



(1.5) 
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where the sum runs over spanning hyperforests F in G, and k{F) is the number of 
connected components of F (note that the second equality in (1.6) uses the following 
Proposition 2.1). If we set wa = 1 for all the hyperedges A & E^-we get as coefficient 
of in the polynomial on the right hand side of the previous equation the number 
of unrooted hyperforests of the hypergraph G with p components. 
If, on the other hand, we specialize (1.4) to A = 0, we obtain: 



/ 



'Dv,t{ip,ip) exp 



E 

AeE 




(1.7) 



where the sum runs over spanning hyperforests F in G with components Fi, . . . ,Fi, 
and V{Fa) is the vertex set of the hypertree F^. This is the generating function of 
rooted spanning hyperforests, with a weight wa for each hyperedge A and a weight ti 
for each root i. By taking the derivatives with respect to ij^, • • • , at t = we easily 
get the generating function of spanning hyperforests rooted at the vertices ii, - ■ • ,ir, 
which is 



/ 



X>y(V','0) (V'V')n • • • (V'V')v exp 



J^WAf, 



(0) 
A 



AeE 



(1.8) 



where we used the shortened notation := ipiipi- If we now set wa = 1 for all 

the hyperedges A G we get the number of hyperforests of the hypergraph G with 
connected components r hypertrees rooted at the vertices ii, ■ ■ ■ , ir- Remark that in 
the case of an ordinary graph f^^'' = = {i'i — — i'j) is a quadratic form in 
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the Grassmann fields, and the previous integral reduces to the evaluation of a reduced 
determinant of the Laplacian matrix, in agreement with the matrix-tree theorem. 

In what follows we shall obtain exphcit formulas for the case of the hypergraph 
fCn which is complete in hyperedges of all possible cardinality, with weight Wg on 
the hyperedges of cardinality s for s — 2,...,n, that is with wa = w\a\- These 
results could in principle and in many cases have been already derived by using the 
standard methods of enumerative combinatorics, that is Lagrange inversion formula 
in connection with the formalism of the exponential generating functions. We hope 
to convince the reader that also in these cases our Grassmann formalism provides 
an alternative, simple and compact way to recover the total weights for rooted and 
unrooted hyperforests on n labeled vertices, which is to say spanning on the complete 
hypergraph /C„. 

This paper is organized as follows. In Section 2 we recall relevant notions from 
graph theory. In Section 3 we illustrate how, at least in the case of complete graph, 
the representation for the generating function of unrooted hyperforests (1.6) can be 
deduced from that for the rooted hypcrtrccs (1.7). In Section 4 we collect all the 
explicit Grassmann integrals that will be used in the following. In Section 5 we show 
the relation between our Grassmann integrals and the explicit solutions achieved by 
standard methods. In Section 6 we deal with rooted hyperforests, while Section 7 
is devoted to unrooted hyperforests. By restricting our general model to the case in 
which only one weight in nonzero, that is Wp = Sp^k , we obtain the explicit evalua- 
tion of the number of rooted and unrooted spanning hyperforests on the /c-uniform 
complete hyper graphs )Cn^ with n- vertices. These results are presented respectively 
in Section 6.1 and Section 7.1. Here we also derive a novel general simple expression 
for the number of unrooted hyperforests with p hypertrees in terms of associated 
Laguerre polynomials and its asymptotic expansion for large number of vertices. We 
consider also another special case, the one in which all the weights are equal, that is 
Wp — 1 for all p, in Section 6.2. We give in Section 7.2 the evaluation of the number 
of hyperforests rooted on p vertices for the hypergraph /C„. Some conclusions are 
presented in Section 8. 

Appendix A collects some basic features of Stirling numbers of the second kind. 
Bell numbers and Bell polynomials. In Appendix B we report, for reader conve- 
nience, the derivation of the number of (hyper-)trees in a unified way by the standard 
exponential generating function formalism and Lagrange inversion formula. In Ap- 
pendix C we provide some results on the asymptotic behaviour of the associated 
Laguerre polynomials which are used in the main text. 

2 Graphs and hypergraphs 

A (simple undirected finite) graph is a pair G = {V,E), where ^ is a finite set 
and is a collection (possibly empty) of 2-element subsets of V. The elements of 
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V are the vertices of the graph G, and the elements of E are the edges. Usually, in 
a picture of a graph, vertices are drawn as dots and edges as lines (or arcs). Please 
note that, in the present definition, loops ( O ) and multiple edges ( <Z>) are not 
allowed. We write \V\ (resp. |£'|) for the cardinality of the vertex (resp. edge) set; 
more generally, we write IS"! for the cardinality of any finite set S. 

A graph G' = {V, E') is said to be a subgraph of G (written G' C G) in case 
V' Q V and E' C E. If V = V, the subgraph is said to be spanning. We can, by a 
slight abuse of language, identify a spanning subgraph {V,E') with its edge set E'. 

A walk (of length A; > 0) connecting vo with Vk in G is a sequence 

(VO, 61,^1, 62,^2, • • ■,ek,Vk) 

such that all Vi e V, all e E, and t'j-i, Vi & Ci for 1 < i < k. A cycle in G is a walk 
in which 

(a) Vq, . . . , Vk-i are distinct vertices of G, and Vk — Vq 

(b) Ci, . . . , Cfc are distinct edges of G; and 

(c) k > 2.1 

The graph G is said to be connected if every pair of vertices in G can be connected 
by a walk. The connected components of G are the maximal connected subgraphs of G. 
It is not hard to sec that the property of being connected by a walk is an equivalence 
relation on V, and that the equivalence classes for this relation are nothing other 
than the vertex sets of the connected components of G. Furthermore, the connected 
components of G are the induced subgraphs of G on these vertex sets.^ We denote 
by c{G) the number of connected components of G. Thus, c{G) = 1 if and only if G 
is connected. 

A forest is a graph that contains no cycles. A tree is a connected forest. (Thus, 
the connected components of a forest are trees.) It is easy to prove, by induction on 
the number of edges, that 

\E\ - \V\ + c{G) > (2.1) 

for all graphs, with equality if and only if G is a forest. 

In a graph G, a spanning forest (resp. spanning tree) is simply a spanning subgraph 
that is a forest (resp. a tree). We denote by J^{G) [resp. T(G)] the set of spanning 
forests (resp. spanning trees) in G. As mentioned earlier, we will frequently identify 
a spanning forest or tree with its edge set. 

^ Actually, in a graph as we have defined it, all cycles have length > 3 (because ei ^ 62 and 
multiple edges are not allowed). We have presented the definition in this way with an eye to the 
corresponding definition for hypergraphs (see below), in which cycles of length 2 are possible. 

^ If V^' C V, the induced subgraph of G on V' , denoted G'[y], is defined to be the graph [V' ,E') 
where E' is the set of all the edges e £ E that satisfy e CV' (i.e., whose endpoints are in V). 
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A rooted tree is a tree with a distinguished vertex called the root. A rooted forest 
is a graph whose connected components are rooted trees. 

Hypergraphs are the generalization of graphs in which edges are allowed to con- 
tain more than two vertices. Unfortunately, the terminology for hypergraphs varies 
substantially from author to author, so it is important to be precise about our own 
usage. For us, a hypergraph is a pair G = {V, E), where V is a finite set and E is a. 
collection (possibly empty) of subsets of V, each of cardinality > 2. The elements of 
V are the vertices of the hypergraph G, and the elements of E are the hyperedges (the 
prefix "hyper" can be omitted for brevity) . Note that we forbid hyperedges of or 1 
vertices (some other authors allow these). ^ We shall say that A ^ E is a k-hyperedge 
if A is a Ar-clcmcnt subset of V. A hypergraph is called k-uniform if all its hyperedges 
are /c-hyperedges. Thus, a 2-uniform hypergraph is nothing other than an ordinary 
graph. 

The definitions of subgraphs, walks, cycles, connected components, trees and 
forests given above for graphs were explicitly chosen in order to immediately gener- 
alize to hypergraphs: it suffices to copy the definitions verbatim, inserting the prefix 
"hyper" as necessary. The analogue of the inequality (2.1) is the following: 

Proposition 2.1 Let G — {V,E) he a hypergraph. Then 

Y.{\A\-1) - \V\+c{G) > 0, (2.2) 
AeE 

with equality if and only if G is a hyperforest. 

Proofs can be found, for instance, in [18, p. 392, Proposition 4] or [17, pp. 278-279, 
Lemma] . 

Please note one important difference between graphs and hypergraphs: every 
connected graph has a spanning tree, but not every connected hypergraph has a 
spanning hypertree. Indeed, it follows from Proposition 2.1 that if G is a fc-uniform 
hypergraph with n vertices, then G can have a spanning hypertree only if /c — 1 divides 
n — 1. Of course, this is merely a necessary condition, not a sufficient one! 

The hypergraph JCn^ has |y| = n vertices and is complete in the /c-hyperedges, in 
the sense that it is /c-uniform and for all choices of k different vertices ii, . . . , in V 
the hyperedge {ii, . . . belongs to the set of its hyperedges. 

The hypergraph^ /C„ has |V^| = n vertices and is complete in the /c- hyperedges for 

all 2 < A; < n, so that E = UL2 E (/ci^^) . 

^ Our definition of liypcrgrapli is the same as tliat of Mc;Cammond and Meier [16]. It is also the 
same as that of Gessel and Kalikow [17], except that they allow multiple edges and we do not: for 
them, E is a, multiset of subsets of V (allowing repetitions), while for us is a set of subsets of V 
(forbidding repetitions). 

*We don't use the symbol /C„ because this is usually used for the complete graph ICn^ with n 
vertices. 
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3 Exponential generating function for hypertrees 
and hyperforests 

Let us consider the complete hypergraph /C„ for every n, with general hyperedge- 
weights wa which vary only with the cardinality of the hyperedge A, i.e. wa — w\a\- 
The /c-uniform complete hypergraph ICn^ corresponds to the case in which the only 
non- vanishing weight is Wk- 

Let tn be the total weight of rooted hypertrees in the case of n vertices \V\ — n, 
w = {wk}k>2 and let 

T{z) = T{z,w):=J2tn{y^)-^ (3-1) 



n>0 



the exponential generating function for the sequence {tn}- 

The exponential generating function for rooted hyperforests is therefore e*^(^\ 
where t counts the number of connected components. 

In the case of complete hypergraphs we can also consider the exponential gener- 
ating function for unrooted trees 



U{z) = U{z,w) :=5^M„(w) 



n>0 



Z" 



(3.2) 



where Un is the weight of unrooted trees in the case of n vertices |V| = n. Of course 
as the root of a trees on n vertices can be chosen in n ways 



and therefore 



and conversely 



Jo ^ 



(3.3) 
(3.4) 
(3.5) 



By using a recursion relation, as it is done in Appendix B, we see that the exponential 
generating function for rooted hypertrees satisfies the relation 



T{z) — z exp 



E 

.k>2 



Wk 



T{z) 



k-l 



{k-l)\ 



(3.6) 



therefore 



z — T exp 



k>2 



{k-l)\ 



(3.7) 
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and by changing variables from lu to T{u!) in the integral in (3.5) we easily get 



U{z)=T{z) + J2^,{l-k)^ 

k>2 



{31 



that is the exponential generating function for unrooted hypcrtrees can be expressed 
in terms of the exponential generating function of rooted hypertrees [3] . Furthermore, 
the exponential generating function for unrooted hyperforests, with parameter A to 
count the number of connected components, is simply e^^^^^ and this can also be 
expressed in terms of the exponential generating function of rooted hypertrees by 
means of (3.8). 

Let us use these relations in order to re-obtain, at least in the case considered here 
of the complete hypergraph, the generating function of unrooted hyperforests in the 
Grassmann representation from the generating function of rooted hyperforests. 

Formula (1.7) for the generating function of unrooted hyperforests for /C„, Siiti — t 
for every vertex, means that 



n\ [z^] e*^(^) = y"p„(V^,V^)exp 
where we shortly denoted 

and 



(0) 
A 



n „ 

n/ 

i=l 



A^E 



(3.9) 



(3.10) 



(3.11) 



iev(!Cn) 

And, if f{z) can be expanded in powers of z, [z"] f{z) is the coefficient of z in the 
expansion. 

It follows that ,for every power r, the coefficient of f is equal to 



n! [z""] T{z) 



T^n{^,i^) {ip,'4^y exp 



(0) 
A 



LA6-E 



(3.12) 



and therefore for each function L defined by a formal power series 



[z^] L[T{z)] = j P„(V^,Vi)Lp,V')] exp 



A&E 



(3.13) 
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Now, the exponential generating function for unrooted hyperforests is e^^^^\ 
where A counts the hypertrees in the hyperforests and we know by (1.6) that 

n\ [z-] e^^(^) = / explx{i^,ij) + J2w\A\ff\ (3.14) 

I AeE ) 



but 



AeE 



(3.15) 



and 



j2^\A\{i-\A\)TA = J2^>'i^-k) E ^^ = E^'^(i-^)%^ (3-16) 

AeE k>2 A:\A\=k k>2 



AeE 

SO that 



J T>n{^,^p) exp i A 



(^,^) + E«^fc(i-^) 



k>2 



+ E^h/^7 (3-17) 



AeE 



But this is exactly formula (3.13) when 

L{y) = e^^(^) 

with 

X(y) :=7/ + E«^fe(l-^) 

fc>2 

which is such that U{z) = K[T{z)\ by (3.8). 



A;! 



(3.18) 
(3.19) 



4 Useful Lemmas on Grassmann integrals 

In the following we shall make use of very simple results for Grassmann integrals. 

Lemma 4.1 Let \ V\ = n be the number of vertices, then 



I 



s\ 



Proof. It trivially follows from induction in n. □ 
We soon derive, by expansion in powers, that 
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Corollary 4.2 Let g be a generic function of the scalar product, that is a polynomial 
as the scalar product is nilpotent of degree n, then 

J V^iij, ij) g {{ij, ^)) = n! g{z) = ^ ^ g{z) 

where the contour integral is performed in the complex plain constrained to encircle 
the origin. 

These are the ingredients to observe that 
Lemma 4.3 Let \V\ = n be the number of vertices, g a generic function, then 

J ^) i^p^p),, ■ ■ ■ (^^),, g ((^, ^)) = ^^^^ j ^) 'tPY g ((^, ^)) 

= (n-r)! [z^^-^] g{z) 

Proof. By integrating over ipi-^jipi-^, ■ ■ ■ on the left hand side we get an 

integral of the form used in the previous Lemma, where both the integration measure 
and the scalar product were restricted on the remaining n — r vertices, so that 

I ^) g {(^, V')) = (n - r)! [z^-^] g{z) . 

By expanding instead on the right hand side we get 

^(^^_^ / Vr,{^,i,){i,,i:r+^[z^]g{z)= (n-r)! [z^-^-] g{z) 

s>Q 

and we get our result by using the previous Lemma 4.1. □ 
Let J the matrix with unit entries for each i,j e V 

Jij^l. (4.1) 

Our common tool is the following 

Lemma 4.4 Let \ V\ =n be the number of vertices, g and h generic function, then 
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Proof. Let us expand the second part of the exponential 

s 

V^) ^)'^ eKW.'^)) [1 + (v;, ^) ^ ((^, ^))] 



^/ 

because all higher powers of (■0, JV') vanish. We get the final hne because in the rest 

of the integral for each i the field -ipi is always multiplied by the companion ipi and thus 
the only contribution in Jif)) comes from the diagonal part, that is □ 



5 Relation with previous approaches 

In virtue of our Lemmas, the Grassmann integrals for the generating functions 
of rooted and unrooted hyperforests at fixed number of vertices can be expressed 
as a unique contour integral of a complex variable. In this section we will show 
the change of variables which explicitly maps those integrals into the coefficient of 
the corresponding exponential generating function in the number of vertices, without 
using the Lagrange inversion formula. 

The sum on all the edges appears in both main formulas (3.9) and (3.17) and in 
our model it becomes 



k>2 A:\A\=k 



(0) 
A 



k>2 
k>2 



{n-k + 1)— — - (V-, (J - Im 



{k-l)\ 



{k-2)\ 



{k-l)\ 



{k-2)\ 



and according to Lemma 4.4, for any function h of the scalar product (V^, -0) 



(5.1) 



/ ^n(0,'0)^ ((0,-0)) exp Yw\A\f^ 

J r>„(0,'0)/i ((0,0)) exp 



(0) 



n^Wk 

k>2 



{k-iy. 



Wk 



k>2 



(^,^)^-^ 

(^-2)! 



(5.2) 
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so the Grassmann integrals reduces to what has been formally obtained in Corol- 
lary 4.2 and we have for (3.9) 



1-E 



k-l 



k>2 



{k-2)\ 



X exp 



,fe-i 



n! / 



2ni J ^"+1 
which is nothing but 



k>2 



ni tT{z) 



k>2 
k-l 



{k-2)\ 



1 



exp 



dz 



(k-iy. 

t^ + ny^/Wk 



k>2 



,tT(z) 



2 TT i / 2;"+! 

with the change of variables (3.6) with T{z) — ^, as 



(5.3) 



(5.4) 



dz dC^ 



Wk 



ifc>2 



(^-2)! 



(5.5) 



Analogously for (3.17) 



1-E 



Wk 



k-l 



k>2 



{k-2)\ 



X exp 
n! / 



A ((^,v^) + 5:^.(1-^)%-^' 

fe>2 



k>2 



{k-i)\ 



2niJ e 



in+l 



e 



fc-i 



k>2 



{k-2)\ 



X exp 



fc>2 



k-l 



k>2 



{k-l)\ 



(5.6) 



which, by using the same change of variables, is nothing but 



\z-\e 



dz 



2 7ii J ^'^+1 
dz 
2 7ri T 



,XU{z) 



exp < A 



T{z)+J2^k{l-k)'^ 

k>2 



(5.7) 
(5.8) 
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6 Rooted hyperforests 



Let us begin with the evaluation of (1.8), that is the case A = which evaluates 
the weight of rooted hyperforests on r vertices ii, ■ ■ ■ ,ir, which on the complete 
hypergraph /C„ does not depend on the particular choice of the vertices, and we 
denote this weight by Vnr- We have 



exp 
{n — r)\ 



k>2 



ni 



exp 



Wk- 



k-1 



k>2 



(k-2)\ 



k>2 



{k-1) 



Of course there are (") different choiches for r different vertices, therefore, if we denote 

by 

En{t;w) =n\[z'']e'^^'^ (6.1) 

the generating function of rooted hyperforests on n vertices, its Grassmann represen- 
tation is 



1-E 



Wk- 



k-1 



k>2 



{k-2)\ 



The expansion in power series of t 



(6.2) 



(6.3) 



r>0 



provides the total weight of rooted hyperforests with r connected components 

tn,r = in,.(w) = ul [z^] [f] 6*^^^) = [f] E^{t- w) (6.4) 

then 



n 



n,r = / T^n{i>,i>) 



ri 



Wk- 



k-1 



k>2 



(k-2)\ 



while then the total weight of rooted hyperforests 

En{w) := K(l; w) = ^t„,^(w) 



|-z^fc>2«'fc (k-iy. (6.5) 



(6.6) 



r>0 
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is given by the generating function at t = 1. 
Let us now introduce the function 



9{x, y; w) := exp 



,fc-i 



k>2 



y 



:=^P,(a;;w) 



s>0 



r 



(6.7) 



which is the exponential generating function for the exponentials Ps{x]w) in the 
variable x, which varies with the choice of the weights w. We recognize that 



X dy 



k>2 ^ ' s>\ 

Therefore the integral in (6.2) can be re-expressed by using 



s- 1)! 



(6.8) 



y 



r+s 



s>0 r>0 



r! S\ 



and 



fe>2 



s>0 r->0 



r! s! X 



(6.9) 



(6.10) 



The same expression could be written also with the help of the derivative with respect 
to the variable let D = 4, then 



E 

fe>2 



,jfc-l 



{k-2)\ 



e'yeix,y;w) 



E 

k>2 



Wk 



(k-2)\ 



D''-^e^ye{x,y;w) 



y 



r+s 



k>2 



{k-2)\ 



s>0 r>0 
r+s 



= En(-;w)E VE 



r! s! 

Wfc 1 



s>0 



{k-2)l[r-{k-l)]l 



E E ^.-M(^;w)J]r 



r! [s- (A; - 1)]! (fc-2)! 



(6.11) 
(6.12) 

(6.13) 
(6.15) 



so that by comparing term by term in (6.10) and (6.15) we recover a recursion relation 
for the polynomials Ps{x, w) 



k>2 ^ ^ 



(6.16) 
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In terms of the polynomials Ps{x,w) we soon get for the generating function of 
rooted hyperforests 



s>0 r>0 



r! si 



1 - 



= E 



r>l 



n — 1 
r — 1 



Pn-r{n;w)f 



(6.17) 
(6.18) 



Therefore the total weight of rooted hyperforests is 



r>l ^ ^ 



(6.19) 



and the total weight of rooted hyperforests with r hypertrees is 



n — 1 
r-1 



Pn-r{n;w) 



from which in particular we obtain for r = 

tn,o — 



(6.20) 



(6.21) 



for all choices of the weights w, a generalization of what occurs for the case of ordi- 
nary trees because the determinant of the weighted Laplacian on the graph is always 
vanishing. 

Also, as Po{x; w) = 1 for all choices of the weights w, of course 



tn.n — 1 



(6.22) 



as there is only one possible hyperforest with n hypertrees, the trivial one in which 
each hypertree is a vertex. 

The weight of rooted hypertrees tn is given by the case r — 1 



tn := tn,i = P„-i(n; w) 



(6.23) 



A more explicit expression for the polynomials Ps{x; w) is obtained by expanding 
the exponential in the definition (6.7) 



Ps{x;w) = s\[y']9{x,y;w) 



j>2 L 



XWi 



y 



ii(i-i) 



n 

.j>2 



XW-, 



W- V(j-i)! 
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so that if we define tlie coefficients Ps,i ( w) by 

P,{x;w) = J2ps,iM^' (6-24) 

1>Q 

we get 

Ps,i = Ps,i{^) = ^' [y'l [^1 ^(^^ w) 

In order to understand the constraint which is imposed in the sum on the coefficients 
/j's, remember that from Proposition 2.1, if Ij is the number of hyperedges of cardi- 
nahty j, n is the number of vertices and r is the number of connected components, 
which in our case is the number of hypertrees 

= - 1) - \V\ + c{G) = J2 W - 1) - n + r (6.25) 

A&E j>2 

and this is exactly the constraint which is imposed. The number / is instead Ue the 
total number of hyperedges. 



n/,!((i-i)!) 



6.1 On the /c-uniform complete hypergraph 

;te hypergraph ICn ^ the hyperedg 

k. This is therefore the particular case of our model in 

(6.26) 



In the /c-uniform complete hypergraph ICn ^ the hyperedges are all the subsets 
A G V of k vertices: \A\ 



which if we introduce the vectors such that their components are 
we have weights 



w — wek 



(6.27) 



and as we wish to count configurations we have to set w = 1 so that in the general 
formulas Wk — 1 and all the others weights for the hyperedges have to be set to zero. 
We have 



9{x, y; Bk) = exp 



X 



y 



k-l 



{k-l)\ 



(6.28) 



and therefore 




x^-i if s = l(k — 1) for integer I 
otherwise 



(6.29) 
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which satisfy the recursion relation (6.16) which for w = takes the form 

Ps{x;ek) = ^ (^^ _ 2) Ps-{k-i){x;ek) . (6.30) 

We easily get that 

P.z(e,) = |w if « = i(^-l) for integer ; ^^^^^ 
I otherwise 

On ICn \ the numbers Ue = I of hyperedges and the number of connected components 
c{G) = r are related by (6.25) 

l{k-l)-n + r^O (6.32) 

that is 

n — r 

ue^I = ^— Y (6.33) 

is the number of hyperedges (of degree k). 

For the number of rooted hyperforests with r hypertrees on the /c-uniform complete 
hypergraph }Cn^ , we have when n — r can be divided by A; — 1 

tnA^k) = Pn-rin;ek) (6.34) 

— l\ (n — r)\ 



r-1 



(6.35) 



feT)![(A:-l)!]^ 
+ [(^-^)"-]l^p-l)n. + rr (6.36) 



r-1 ynglp-l)! 
where the prefactor in (6.35) 

{n — r)\ 



(f5f)!P-l)!]- 



n — r 
1 



(6.37) 



is exactly the number of ways in which n — r vertices can be divided into {n—r) /{k — 1) 
groups of A; — 1 elements and in (6.36) we have replaced the dependence from the 
number of vertices n with that from the number of hyperedges ue- 
In the case of simple graphs {k — 2) it follows that 

W(e2)= (6.38) 
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which at r = 1 provides the well-known result by Cayley about the number ti„ of 
spanning unrooted trees on the complete graph with n vertices 



n 



Also 



which could be obtained by direct evaluation as 



n-l 



(6.39) 
(6.40) 



1 - 



n 



n + t 



This relation says at t = 1 that the total number of rooted forests is 

K(e2) = (n + 1)"-^ 



(6.41) 



(6.42) 



In this simple case also the whole generating function can be expressed in terms of 
the generalized exponential [1] (the usual exponential is at a = 0) 



^a(^) := ^ (an + 1)" ^ 



n>l 



Z 



which satisfies 

Saiz)'" In Sa{z) = z Sa{z) = S{az)'' 
where E{z) is a shorthand for £i{z). Indeed 



(6.43) 



(6.44) 



n-l 



n>l 



n>l 



^Si{tz) =S{zY = e''^^'^ 
6.2 On the complete hypergraph 



(6.45) 



We shall consider here the complete hypergraph /C„ when all the hyperedge- 
weights Wd are set to one, that is 

w = 1 (6.46) 



where 1 is the vector with 1 on all components. We have 

0{x,y;l) :=exp[x(e^-l)] = 



s>0 



s\ 



(6.47) 
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where bs{x) are the Bell polynomials, see Appendix A, and therefore 



l>0 



so that 



(6.48) 



(6.49) 



where | ^ } is a Stirling number of the second kind, and it is the number of ways to 
partition a set of cardinality s into I nonempty subsets. 
The recursion relation (6.16) becomes here 



k>l ^ ^ 



(6.50) 



The number of rooted hyperforests with r hypertrees on the /c-uniform complete 
hypergraph ICn is therefore 



n 



TIE 



nE>0 



n — r 



He 



(6.51) 



K{n) 



(6.52) 



and the total number of rooted hyperforests is 

Enil) = 5^ W(l) = E I j) ^-^W = 
r>l k>l ^ ^ 

because of (6.50) for x — n. 

7 Unrooted hyperforests 

According to our general formula the generating function for unrooted hyperforests 
on n vertices is given by the Grassmann integral 



F„(A;w) :=n! [z^] e^^^'^ = 

Vr,{i!,2p) exp { A 



(^,^) 



X exp 



k>2 



k-1 



k\ 



k>2 



{k 



fc-2 



k>2 



{k-2)\ 



which we expand in A 



F„(A;w) = E^-.p(w)^'' 

p=0 



(7.1) 



(7.2) 
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where 'U„,p(w) is the total weight of unrooted hyperforests with p hypertrees. 

We find convenient to introduce the polynomials ns(A;w) and the coefficients 
7rs,r(w) according to 



exp 



k>2 

It soon follows that 



w 



s>0 



s>0 p>0 



r 

s! 



(7.3) 



F„(A;w) = J]n,(A;w) /"Pn(V', 

S>1 



exp 



fc-2 



{k-1) 



k>2 



= X]n^(-^; w)t„,^(w) 

= E(^ZJ) n.(A;w)P._,(n;w) 



(7.4) 
(7.5) 

(7.6) 



The total weight of unrooted hyperforests on the set on n vertices, irrespective from 
the number of hypertrees, is obtained from the partition function at A = 1 



F,(w) := F„(l; w) = ^ n,(l; w) t„,,(w) 



S>1 



Also we get 



Un,p = 1in,p(w) =E 7rs,p(w) in,s(w) 



S>1 



(s - i) ^n-5(^; W) 



Remark that from the definition 

7i"s,p(w) = when p > s 
so that lis (A; w) is a polynomial of degree s. It is monic because 

7rs,s(w) = 1 . 
And remark also that TTsfiiw) — while 



for s = 1 
(1 — s) otherwise. 



(7.7) 

(7.8) 
(7.9) 

(7.10) 
(7.11) 

(7.12) 
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Accordingly 'U„,o(w) = and 'U„,n(w) = 1, while it follows that the weight of unrooted 
hypertrees on n vertices is simply the weigth of the rooted hypertrees divided by n, 
indeed from (7.8) 



,(w) := 'Un,i(w) = Pn_i(n; w) + ^ (1 - s) {^^ _ Pn-s{n\ w) (7.13) 
= P„_i(n;w)- (n-1) Y.'^s {^^Z^^^ Pn-s{n-^) (7.14) 



Pn_i(n;w) 



n 



^n(w) 



n 



(7.15) 
(7.16) 



where we used the recursion relation (6.16) for the polynomials Ps{x]\/) &t x = n 
and s = n — 1. 



More formally we can follow a different strategy. Let D — ^ then 



exp 



k>2 



r 
k\ 



exp 



y^^fc (1 - fc) 



k>2 



exp (t y) 



t=x 



so that 

and therefore 
Pn(A;w) = 

= exp 



ns(A, w) = exp 



A y^Wfc (1 - A;) 



fc>2 



IT 



(7.17) 



(7.18) 



fc>2 



k>2 



fe-i 



k>2 



{k-2)\ 



J 



(7.19) 



: exp 



Jjk ^ D^~^ 

A J2'^kil-k)— + nJ2wkjj^^ 



k>2 



k>2 



k>2 



n — 1 
k-2 



n-k+l 
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now, we expand first the second exponential, to get once more 

Wk (1 - k] 

k>2 



F„(A; w) = exp 



= exp 



A y^Wfc (1 - A;) 



E 

s>0 



k>2 

n — 1 
s-1 



k\ 

r 

ns(A,w) Pn-s{n,w) . 



(7.20) 



S-1 



t=x 



7.1 On the A;-uniform complete hypergraph 

When w = ejfc the formula (7.3) becomes 



exp 



.k\^ 



Ah/+ (1-^) 



k\ 



E n.(A; 



s>0 



5^ ^sAek) X' ^ . (7.21 



s>0 p>0 



We introduce a family of generalized Hermite polynomials ifi*^^ (x) as defined by the 
generating function 



exp 



X z + {1 — k) 



k\ 



s>0 



s\ 



(7.22) 



which when k — 2 are related to the ordinary Hermite polynomials Hg by 



Hi'\x) = lie, 



.X 



22 V 22 



(7.23) 



where He^ are sometimes used [22]. Similar generalizations of the Hermite polynomials 
can be found in [23-25] . We then get 



H,(A;e,) = Ati/f (a*^) . 
Thus the generating function of unrooted hyperforests is 



(7.24) 



Fn{X;ek)= V f"-;) ("-^)' At fff (a^) (7.25) 



p: (n-p)l(fc-l) 
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where the sum is restricted to the values of p such that n — p can be divided by A; — 1. 
By using (7.19) we get instead 



-P'n(A;efc) =exp 

—Xk exp 
\ - 

—Ah exp 



r - n 



n — 1 



n—k+l 



k-2 

n fn — 1 
k-2 



fc-i 
A fc 



j/n—k+1 



A'^ (^-1)! 
In the particular case /c = 2 we soon get 



Kit) - -3T 



n fn — 1 



A ~ 2 



-f^^_fe+i(t) 



F„(A;e2)= Va" 



(• 



(7.26) 



because exp [a ^] is the translation operator from ^ to ^ + a. The same result can 
be obtained by using (7.20) and (6.41) as 



-^n(A;e2) = exp 
= exp 

= exp 

= Va" 



-A 
-A 



2 
"2" 



En{t,e2) 



t=x 



[(t + n)"-n(t + n)"-^] 



t=\ 



t + 



n 



n 



n 



n-l 



t=Vx 



He, 



n 

71 



This formula has been reported in [26] for A = 1, where it counts the total number 
of unrooted forests. In this case (7.25) becomes instead 



^"(^2) = E _ ^"-^He,(l) 



(7.27) 



in agreement with what obtained in [26] and reported as the series A001858 in the 
The On-Line Encyclopedia of Integer Sequences by Sloane [27]. 
By using D = we get 



exp 



X z + {1 — k) 



k\ 







= exp 





exp [x z\ 



(7.28) 
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and therefore 



= exp 

= E 

q>0 

= E 

q>0 



1-k 



k\ 

1 fl-k^'' 



X 



1 n-kv 



X 



si 



{s — kq)\ 



x" 



-kq 



1 fl-k 



{s-kq)\ 



y-{k-l)q 



which imphes because of (7.24) 

n,(A;e,) = J] V: 

so that 

and therefore, by using (7.8) 

^ , \p + q{k-l)]\ 1 [l-k\ 

UnA^k) - 2^tn,p+qik-l)[ek) (^3^)1 ^ J 



I fl-k 



(s — kq)\ 



(5. 



'p,s — {k—l)q 



g>0 

(n- 1)! 



n 



(fc-l)! 



f p\p + {k — l)q f 1 — k 



kn 



(7.29) 
(7.30) 

(7.31) 

(7.32) 
(7.33) 

(7.34) 
(7.35) 



when n — p can be divided by k — 1, otherwise it vanishes, where we used the relation 
(6.20) and the exphcit expression (6.35). Once more in the simpler case k = 2 this 
formula reduces to 

-...(e^) = ^E {-IJ (:) (pi;! i) n'-'-'i. + .y. (7.36) 



a result which can be found in [11,28]. 
In order to proceed we need the sums 



E 



p 



^A<l){v-q)\ dz v\ P 



p 



L(;-^\z) + zL^P 



(v-p+l). 



(7.37) 

(7.38) 

(7.39) 
(7.40) 
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where Lm\x) are the associated Laguerre polynomials 



m + q;\ {—xY 



i/=0 ^ / 



m — 1/ J v\ 

i/=U ^ 

which satisfy the recursion relation 

L%{z)^ \ [pL^^\z)-{p + k)L^\{z) 
We arrive at the representation 
(n-1)! n 



(7.41) 



(7.42) 



{k-l)\ 



71 — p 



k-1 



kn 

kn \ , ^ ^ V t 



kn 
k-1 



(7.43) 



for the number of unrooted hyperforests with p hypertrees on the /c-uniform complete 
hypergraph K.n'' with n vertices. 

In order to study the asymptotic behaviour of the previous expression in the limit 
of large n at fixed p we need the following expansion for the Laguerre polynomial 

that can be easily obtained from the definition (7.41), as shown in Appendix C, then 



pLp 



(S-p) f kn 



k-1 



(S-p) f k 



-n)P 1 k 



k-1 J {p-l)\nk-l 



+ [n - p) Lit-,' 

because the leading terms in the two contributions cancel out. We get 

(n-l\ {n-p)\ n^-^ /fc-1^^"' 

Remark that when p—1 this formula is exact, indeed 



(7.45) 



(7.46) 



\n — 1)' n^- 



(7.47) 



is the number of unrooted hypertrees in n vertices, because of the general result (7.16) 
and the explicit expression (6.35). In [19] this number is quoted as obtained in [20]. 
The formula (7.46) at A; = 2 provides the result 



n 



1\ n'^-P-^ 



p-lj 2P-1 



(7.48) 
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already obtained in [11] by a different method. It follows that the partition function 
is, if A is such that the relevant contribution to the sum comes from regions which 
don't change with n, a problem which we will discuss elsewhere, we get 

|«..(e.)A-«»-A|("-)(A)^„.-..(x,i_)"-' 

which at A = 1 provides the well-known result by [9, 10]. 

More generally, by using the Stirling approximation for large factorials 

UnA^k) ^ —i^ ^7 7T7 —j— 7.49 

while 

Y.-nA^^)>^^^'^ ^^f^. Ae^[(^-^)'l^^ (7.50) 

7.2 On the complete hypergraph 

When w = 1 (7.3) becomes 

exp [A (1 - y) {ey - 1)] = J] n,(A; 1) ^ = E E "^^M !t (^-^l) 

Now 



s>0 ' s>0 p>0 



TTsA^) = 5![|/^][A^']exp[A(l-y)(e^-l)] (7.52) 

(ey - lY 

= s! [/] (1 - Z/)^ ^^p^ (7.53) 



so that the number of unrooted hyperforests with p hypertrees obtained by formula 
(7.8), by using the number of rooted hyperforests given in (6.51), is 

^n,p(l) = (7.56) 

p + g 1 s\ 

'^-^^ ' -n) 1 



i:(::i)^"-.("'|:(-)'-"C-:-, 
E(::J)E{"r}"'5:(-)-(,_:_ 

s>l ^ ^ r>0 ^ ■' q>0 ^ ^ 



p j {p + qV- 

p + q \ s\ 



qJ { p } {p + qV- 



26 



Of course, because of the general result (7.16), 



M„(l) = M„,i(l) = ^ = ^ } ^ < }n (7.57) 

n n ^-^ I r I 

r>0 ^ ^ 

a sequence which is reported with the number A030019 in the The On-Line Encyclo- 
pedia of Integer Sequences by Sloane [27]. 

8 Conclusions 

We have studied the generating function of both rooted and unrooted hyperforests 
in the complete hypergraph with n vertices, when the weight of each hyperedge de- 
pends only on its cardinality. All the results could also be obtained by starting from 
recursion relations in the number of vertices, to obtain implicit relations for the for- 
mal power series of the generating function, which can be afterwards solved by using 
the Lagrange inversion formula. However we showed here how the same problem can 
be directly and more easily solved by means of a novel Grassmann representation. 

Once we obtained the general solutions we have restricted to particular cases to 
recover more explicit results. In particular we considered the case of the /c-uniform 
complete hypergraph, where only edges of cardinality k are present. When this weight 
is set to one we are reduced to a counting problem. We thus obtained a generalization 
of many known results in the case k = 2 namely of ordinary forests on the complete 
graph. In the case of unrooted hyperforests we also recovered a novel explicit expres- 
sion for their number with p connected components, that is hypertrees, in terms of the 
associated Laguerre polynomials, for any k. We have also presented the asymptotic 
behaviour of these numbers for large number of vertices. 

A second direct application of the general solutions is obtained for the complete 
hypergraph when all the hyperedges have the same weights. 
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A Stirling and Bell numbers, Bell polynomials 

The Stirling numbers of the second kind, denoted by | ^ } according to the notation 
introduced in 1935 Jovan Karamata and promoted later by Donald Knuth, stands for 
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the number of ways to partition a set of cardinality n into k nonempty subsets. Thus 

1^1 = ior n< k 

- for n > 1 



1 



if we agree that there is one way to partition an empty set into zero nonempty parts. 
Chosen an object among n > to be partitioned into k nonempty parts, we either 
put it into a class by itself (in { ^l} } ways) or we put it together with some nonempty 
subset of the other n — 1 objects (there are A; {"~^} possibilities, because each of the 
{ } ways to partition the n — 1 other objects into k nonempty parts gives k subset 
that the chosen object can join), hence we get the recurrence 

n — ll (n — 1 



which enables us to compute them. 

Their exponential generating function is 



n>0 n>k 



The Bell number 6„ is the number of all possible subsets of a set of cardinality n, 
hence 

A;>0 

Their exponential generating function is 

n>0 n>0 fe>0 fe>0 n>0 fc>0 

The Bell polynomials, also called exponential polynomials, are given by 

jk>0 

SO that 

hn{l) = &n (A.6) 

and they satisfy the recurrence relation 

hn+i{x)=x[hn{x) + h'^{x)] (A.7) 
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as follows from (A.l). 

Their exponential generating function is 

n>0 n>0 fc>0 ifc>0 

B Exponential generating function 

(2) 

Counting the number of unrooted trees Un on the complete graph ICn is presented 
in [1, Chapter 7] as a simple application of the formalism of the exponential generating 
function. 

For n > the recurrence 

M„=V^ Y] / nl \ a^...a^Ua,---Ua^ (B.l) 

m>0 ai,a2,--- ,am ^ ^ m/ 

aiH ham=n-l 

can be obtained as follows. A given vertex is attached to m components of sizes 
There are " ) ways to assign n — 1 vertices to those components 
and ai - ■ ■ ttm ways to connect the given vertex to them. There are ■ ■ ■ ways 
to connect those individual components with spanning trees; and we divide by m! 
because the m components are not ordered. 
As the number of rooted trees is 

tn = nun (B.2) 
the recurrence relation can be re-written as 

^ = EA E (B.3) 
n! ^ m! ^ ai! a^! 

aiH \-am=n-l 

By introducing the exponential generating function for the sequence 



ni 

n>0 



it follows that the inner sum in (B.3) is the coefficient of in T{z) 



[-"-^] E^^(-r= [z--^]e-^^^ (B.5) 



nl '- ^-^ ml 

m>0 



where we have included also the case n = 1 by adding the contribution m — 0. And 
therefore 

T{z) = ze'^^'^ (B.6) 
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so that because of (6.43) and (6.44) T is related to the generaUzed exponential £ by 

T{z)^z£{z). (B.7) 
T is also related to the Lambert W function [29] by 

T{z) = -W{-z) . (B.8) 

Now 

n n 

This result is usually attributed to Cayley in 1889 [30], but in his paper he refers to 
a previous result by Borchardt in 1860 [31]. 

More generally when 6{u) is a formal power series in u with ^(0) = 1, a relation 
for the formal power series T[z) of the form 

T{z) ^ ze{T{z)) (B.IO) 

has a unique solution, which is given by Lagrange inversion formula [2] 

[^"] T(z) = - [T"-!] ^ (T)" . (B.ll) 
n 

Furthermore 

[2«] T{zy = ^ [T"-"] 9 {T f . (B.12) 
In our application to the trees, 0{T) — and therefore 

n n n n 

While the number of rooted forests with r trees is given by 

tn,r = ^ [^1 TizY = 1^ [T-] = : J) n'^-'- . (B.14) 

(k) 

More generally, in the case of the fc-uniform complete hypergraph ICn with weights 
Wk the recurrence relation for the weight of unrooted hypertrees is 

m 

u ^ T V ( 

m>0 \k-\)- r\ ai,a2,--,am ^ '-^ ' ^ 
m/(fe-l) aiH ham=n-l 

(B.15) 

where at variance with respect to (B.l) the sum on m is restricted to integers that can 
be divided by A; — 1 and appears a combinatorial factor m- because this 

(^)![(fe-l)!]S=T 
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is the number of ways in which the m sub-hypertrees can be hooked to the starting 
vertex by using hyperedges of cardinahty k. As a consequence the equation for the 
rooted hypertrees becomes 



which is to say 



T{z) = ze""^^^ (B.17) 

We can now apply the Lagrange inversion formula with e{T) = e^" (^-D! and there- 
fore 



n-l 
1 



-1 



n n n ' ' n (n^^l P _ 1)1]^ 

(B.18) 

While the weight for the rooted hyperforests with r hypertrees is 



n—r 
I fc-1 



[.Ji(.) [7 Je -(r-1)! (|5f)!p-l)!]f^ 

(B.19) 

when {n — r)/{k — 1) is an integer. It is indeed the total number of hyperedges. 

In the general case of the complete hypergraph /C„ the recurrence relation for the 
total weight of unrooted hypertrees is more involved, but the possibilities of attaching 
hyperedges of different cardinality at the starting vertex are mutually avoiding and 
this makes the recursion affordable. It follows that the generating function satisfies 
the equation 

so that 



(r-1)! 

where we introduced the polynomials P„_r(w) of the weights w^s defined in (6.7). 

In the simpler case in which all the weights are equal to, say, x, the recurrence 
relation for the unrooted hypertrees is 

= ~T IT) ^' ( ^ ^ )ai---amUa^---Ua^ (B.22) 

m>0 l>0 ai,a2,---,am ^ ^' ^ rn/ 



oiH \-am=n-l 



where, at variance with respect to (B.l) there appears a factor |™} because this is 
the number of ways in which the m sub-hypertrees can be hooked to the starting 
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vertex by using I generic hyperedges. As a consequence the equation for the rooted 
hypertrees becomes 

tn r „_n NT^ 1 (m ' 



(B.23) 



m>0 />0 



which is to say 

T{z)^ ze<'^^^'-^) (B.24) 

that is (B.20) for Wk = x for all k, a relation that in the case a; = 1 is reported in the 
Warme's Ph. D. Thesis [21] as due to W. D. Smith, but see also [17]. We can now 

apply the Lagrange inversion formula with 6{T) — e^^'^ and therefore 



n n ^ ' ' n ^ n 

While the total weight of rooted hyperforests with r hypertrees is 



(B.25) 



(n-l)\ 
(r - 1) ! 



^y„-r] ^n.{eT-l) ^ {^^ b^-riu x) . (B.26) 



C Asymptotic behaviour of associated Laguerre 
polynomials 

In this appendix we will study the asymptotic behaviour of the associated Laguerre 

polynomial Ls~^ large n. 

We remark that for a, i^, s all integers and a ^ 1 



{s + ay. 



a'-" + a'-"-^ [s + --- + {u + l)] 



a''" + a 



s—u—l 



s{s + l) u{u + l) 



if 



n — p 
d-1 



-P 



we get, for n ^ 1 at first order in 1/n 

n 



d-1 



l-is-u) 



pd 



n 



and 



[s + ay 
{u + ay 



n 



d-1 



" pd d — 1 
l-{s-u) — + 



n 



n 



s{s + l) u{u + l) 
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so that 



dn 
d-1 



E 

i/=0 



{s + a)\ 1 /_ dn Y 



d 



\ S -J s 



u=0 



V 



, , pd d — 1 



n n 



s(s + l) u{u+l) 



Now 



E (') i-dr = (1 - dy 



and by taking one and two derivatives with respect to — d we get 

i24l)i-dr' - s{i-dr' 



s-2 



and therefore 



s (s - 1) (1 - dy-^ - 2 s (1 - d)^-^ d 

s (s + 1) (1 - dy-^ d^-2s{i- dy-^ d 



and we get 



+ 



d-1 



n 



^(^±1) (1 _ dy - (1 - dy-' d' + s{i- dy-' d 



-n] 



s\ 



1 + 



s{p+l)d s{s + l){l-2d) 



+ 



{d-l)n 2{d-l)n 



from which (7.44) follows. 
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